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Abstract 

In  this  work  we  design  surveillance  trajectories  for  a  network  of  autonomous  cameras  to  detect  intruders  in  an  environment. 
Intruders,  which  appear  at  arbitrary  times  and  locations,  are  classified  as  static  or  dynamic.  While  static  intruders  remain 
stationary,  dynamic  intruders  are  aware  of  the  cameras  configuration  and  move  to  avoid  detection,  if  possible.  As  performance 
criteria  we  consider  the  worst-case  detection  time  of  static  and  dynamic  intruders.  We  model  the  environment  and  the  camera 
network  by  means  of  a  robotic  roadmap.  We  show  that  optimal  cameras  trajectories  against  static  intruders  are  obtained  by 
solving  a  continuous  graph  partitioning  problem.  We  design  centralized  and  distributed  algorithms  to  solve  this  continuous 
graph  partitioning  problem.  Our  centralized  solution  relies  on  tools  from  convex  optimization.  For  the  distributed  case, 
we  consider  three  distinct  cameras  communication  models  and  propose  a  corresponding  algorithm  for  each  of  the  models. 
Regarding  dynamic  intruders,  we  identify  necessary  and  sufficient  conditions  on  the  cameras  locations  to  detect  dynamic 
intruders  in  finite  time.  Additionally,  we  construct  constant-factor  optimal  trajectories  for  the  case  of  ring  and  tree  roadmaps. 
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1  Introduction 

Remote  surveillance  of  human  activities  for  civil  and  military  applications  is  receiving  considerable  attention  from 
the  research  community.  Public  areas  such  as  banks,  art  galleries,  private  houses,  prisons,  department  stores,  and 
parking  lots,  are  now  equipped  with  camera  networks  to  detect  important  activities,  [18,  22].  From  a  technological 
perspective,  one  of  the  main  challenges  consists  of  developing  efficient  algorithms  for  the  cameras  to  autonomously 
and  distributively  complete  tracking,  surveillance,  and  recognition  tasks. 

In  this  work  we  focus  on  the  problem  of  detecting  intruders  by  means  of  a  network  of  autonomous  cameras.  In 
particular,  we  consider  Pan-Tilt-Zoom  (PTZ)  cameras  installed  at  important  locations.  We  assume  the  cameras 
to  move  their  field  of  view  (f.o.v.)  to  cooperatively  surveil  the  whole  environment.  We  develop  algorithms  for  the 
cameras  to  self-organize  and  to  detect  intruders  in  the  environment,  that  appear  at  arbitrary  locations  and  times. 
We  consider  static  intruders,  which  remain  stationary,  and  dynamic  intruders,  which  move  to  avoid  detection,  if 
possible.  As  performance  criteria  we  consider  the  worst-case  detection  time,  that  is  the  longest  time  needed  for  the 
cameras  to  detect  intruders. 

Related  work.  Works  related  to  our  camera  surveillance  problem  can  be  found  in  the  mobile  robotics  and  com¬ 
puter  science  literatures.  In  mobile  robotics,  the  patrolling  problem  consists  of  scheduling  the  motion  of  a  team  of 
autonomous  agents  in  order  to  detect  intruders  or  important  events,  e.g.,  see  [1,  4,  17,  20].  It  should  be  observed 
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that  the  patrolling  problem  and  the  problem  considered  in  this  paper  significantly  differ.  Indeed,  cameras  are  fixed 
at  predetermined  locations,  and  their  f.o.v.s  must  lie  within  the  cameras  visibility  constraints.  On  the  other  hand, 
robots  are  usually  allowed  to  travel  the  whole  environment,  and  are  usually  not  subject  to  visibility  constraints.  Con¬ 
sequently,  algorithms  developed  for  teams  of  robots  are,  in  general,  not  applicable  in  our  setup.  Similarly,  algorithms 
for  graph-clearing  and  graph-search  do  not  extend  to  our  scenario  [13,  14,  19]. 

In  the  context  of  camera  networks,  the  perimeter  patrolling  problem  has  recently  been  studied  in  [4,  8,  21].  In  these 
works,  distributed  algorithms  are  proposed  for  the  cameras  to  partition  a  one-dimensional  environment,  and  to 
synchronize  along  a  trajectory  with  minimum  worst-case  detection  time  of  intruders.  We  improve  the  results  along 
these  directions  by,  for  instance,  developing  cameras  trajectories  and  partitioning  methods  for  general  environment 
topologies. 

In  this  work  we  present  algorithms  for  graph  partitioning.  It  is  worth  noting  that  our  graph  partitioning  problem 
differs  from  classical  setups,  e.g.,  see  [2,  3,  11,  12].  Indeed,  in  these  works  the  graph  partitioning  problem  is  usually 
combinatorial,  and  it  consists  of  partitioning  the  vertices  or  the  edges  as  to  optimize  a  certain  performance  function. 
Instead,  we  formulate  continuous  graph  partitioning  problems,  in  which  the  graph  is  a  physical  entity,  and  the 
partition  is  obtained  by  splitting  the  edges.  As  it  will  be  clear  in  the  sequel,  our  results  on  graph  partitioning 
are  general  and  applicable  to  different  problems.  For  instance,  if  each  edge  of  the  graph  represents  a  task  to  be 
accomplished  by  the  processors  at  its  endpoints,  then  our  algorithms  can  be  used  for  dynamic  load  balancing  for 
multiprocessor  networks  [9,  16]. 

A  preliminary  version  of  this  work  appeared  in  [6].  Extensions  with  this  work  include  (i)  the  design  of  trajectories 
against  static  intruders  for  cyclic  roadmaps,  (ii)  the  design  of  trajectories  against  dynamic  intruders,  and  (iii)  a 
characterization  of  necessary  and  sufficient  conditions  for  finite  detection  time  of  dynamic  intruders. 

Paper  contributions.  The  main  contributions  of  this  work  are  as  follows.  First,  we  propose  the  continuous  graph 
partitioning  problem,  in  which  a  partition  of  a  weighted  graph  is  obtained  by  splitting  the  graph  edges,  and  the  cost 
of  a  partition  equals  the  longest  length  of  its  parts  (Section  2).  We  show  that  the  continuous  graph  partitioning 
problem  is  convex  and  non-differentiable,  and  we  characterize  its  solutions.  Then,  we  derive  an  equivalent  convex 
and  differentiable  partitioning  problem,  which  is  amenable  to  distributed  implementation. 

Second,  we  define  the  camera  surveillance  problem  for  the  detection  of  static  and  dynamic  intruders  (Section  3).  We 
model  the  environment  and  the  camera  network  by  means  of  a  robotic  roadmap,  and  we  formalize  the  worst-case 
detection  time  of  static  and  dynamic  intruders. 

Third,  we  exhaustively  discuss  the  case  of  static  intruders  (Section  4).  We  show  that,  for  tree  and  ring  roadmaps, 
cameras  trajectories  with  minimum  worst-case  detection  time  can  be  designed  by  solving  a  continuous  graph  par¬ 
titioning  problem.  For  general  cyclic  roadmaps,  our  trajectories  based  on  continuous  partitions  are  proved  to  be 
optimal  up  to  a  factor  2.  However,  we  conjecture  that  optimality  is  achieved  also  in  this  case. 

Fourth,  for  the  case  of  dynamic  intruders,  we  derive  a  necessary  and  sufficient  condition  on  the  cameras  locations 
for  the  existence  of  a  trajectory  with  finite  detection  time  (Section  5).  We  focus  on  ring  and  tree  roadmaps.  In 
particular,  for  the  case  of  ring  roadmaps  we  design  a  trajectory  with  detection  time  within  a  factor  3/2  of  optimal. 
Instead,  for  tree  roadmaps,  the  performance  of  our  trajectory  is  within  a  factor  2  of  optimal. 

Fifth  and  finally,  we  consider  three  different  communication  models,  and  we  propose  distributed  algorithms  for 
the  cameras  for  continuous  graph  partitioning  in  all  these  scenarios  (Section  6).  In  particular:  our  first  algorithm 
assumes  a  synchronous  mode  of  operation  of  the  cameras;  our  second  algorithm  assumes  an  asymmetric  broadcast 
communication  model  and  extends  the  class  of  block-coordinate  descent  algorithms  to  the  constrained  case;  and  our 
third  algorithm  only  requires  gossip  communication.  We  prove  convergence  of  all  these  algorithms,  and  we  analyze 
their  performance  in  a  simulation  study. 


2  Continuous  Partitions  of  Weighted  Graphs 

In  this  section  we  introduce  the  problem  of  continuous  graph  partitioning.  A  solution  to  this  problem  will  be  used 
to  design  optimal  cameras  trajectories. 
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Let  G  =  (V,£)  be  an  undirected  weighted  graph,  where  V  and  £  denote  the  vertex  and  edge  sets,  respectively.  Let 
£ij  G  R>o  be  the  weight  associated  with  the  edge  {vi,Vj}  G  £.  For  a  subset  of  vertices  Vc  C  V,  define  the  i-th  set  of 
neighbors  as  AT,  =  A f™  U  Af°ut ,  where 


7V‘n  =  {vj  G  Vc  :  {vi,Vj}  G  £}, 
Af°ut  =  {vj  G  V\VC  :  {Vi,Vj}  €£}. 


A  continuous  partition  of  the  weighted  graph  Q  is  a  set  P  =  {Pi, . . . ,  Vn },  where  (see  Fig.  1) 

Pi  |^J  [Vi,  Vij], 

Vj^Mi 


(2) 


and  Vij  G  [i>i,  Vj]  is  defined  by  some  G  [0, 1]  as  1 


!Vi  +  a,ij  (vj  -v^,  if  j  G  ATim,  i  <  j, 

Vi  +  (1  -  aji){vj  -  Vj),  if  j  G  A 4m,  i  >  j,  (3) 

vj,  ifjGMout. 

The  dimension  of  a  continuous  partition  equals  £(P)  =  max{Li, . . . ,  Ln },  where  Li  is  the  sum  of  the  lengths  of  the 
segments  in  Pi,  i  =  1, . . . ,  n.  In  other  words 

Li—  'y  '  cxijlij  +  ^  ~  (1  ~  aji)£ij  +  £ij- 

Vj  eM-"  ,i<j  Vj  c.Vy  ,i>j  v,  cArf''"1' 

Let  L  and  a  be  the  vectors  of  Li  and  a.ij,  respectively.  Notice  that  a  continuous  partition  is  entirely  specified  by  a 
parameters  vector  a. 

Let  Q(V c)  =  ( Vc,£c )  be  the  subgraph  of  Q  induced  by  the  vertices  Vc,  where  £c  =  (Vc  x  Vc)  <~)£ .  Define  the  weighted 
incidence  matrix  A  G  RlVclxl£d  as 


!lij,  if  e  =  {vi,  Vj}  G  £c,  i  <  j, 

-t ij ,  if  e  =  {vi,Vj}  G  £c,  i  <  j,  (5) 

0,  otherwise, 

and  the  weighted  incidence  vector  b  G  as 

bi  =  y  tij  +  y '  £ij.  (6) 

VjGj\Tla,i>j  Vj^M^ 

Notice  that  L  =  Aa  +  b,  and  that  for  every  a  G  it  holds 

||Aa  +  6||i=  ^  tij. 

{vj,vj}e£ 


Let  0  and  1  be  the  vectors  of  all  zeros  and  ones,  respectively.  We  address  the  following  minimization  problem. 

Problem  1  (Continuous  min-max  partition)  Given  a  weighted  graph  Q  =  (V ,£)  and  a  subset  of  vertices  Vc  C  V, 
let  A  and  b  be  as  in  (5)  and  (6),  respectively.  Determine  a  continuous  partition  a ^  satisfying 

11^0:^+61100=  min_||Ao  +  6||oo,  (7) 

ot<.ot<.ot 


1  Given  any  two  points  Xi,  Xj  G  Rm  for  some  m  G  N,  we  let  [xt,  Xj]  =  {tXi  +  (1  —  t)xj  :  t  €  [0, 1]}. 
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for  some  constraints  vectors  a,  a  with  0  <  a  <  a  <  1. 


Notice  that  (7)  is  a  convex  minimization  problem,  for  which  efficient  centralized  solvers  exist  [7].  On  the  other 
hand,  since  (7)  is  not  differentiable,  distributed  solvers  may  be  difficult  to  implement.  We  next  derive  an  equivalent 
differentiable  minimization  problem,  which  is  instead  amenable  to  distributed  implementation. 

Problem  2  (Continuous  min  partition)  Given  a  weighted  graph  Q  =  (V,  8)  and  a  subset  of  vertices  Vc  C  V,  let 
A  and  b  be  as  in  (5)  and  (6),  respectively.  Determine  a  continuous  partition  at,  satisfying 

|| -Ac*?!  +  foj|  2  =  min_  \\Aa  +  foj|2.  (8) 

at<.at<OL 


for  some  constraints  vectors  a,  a  with  0  <  a<  a  <  1. 

Remark  3  (Uniqueness  of  partitions)  Since  the  minimization  problem  (8)  is  strictly  convex,  the  continuous  min 
partitioning  problem  admits  a  unique  minimum  value,  and  the  set  of  minimizers  is  a  singleton  if  and  only  if  the 
matrix  A  has  a  trivial  null  space.  It  can  be  shown  that  A  has  a  trivial  null  space  if  and  only  if  the  induced  graph  Q(VC) 
is  a  tree.  In  particular,  if  the  graph  G(VC)  is  connected,  then  the  dimension  of  the  null  space  of  A  equals  \8C\  —  |  Vc|  + 1. 

Remark  4  (Unconstrained  partitions)  Consider  the  (unconstrained)  partitioning  problems  (7)  and  (8)  with 
a,  a  £  RI£cI.  It  can  be  verified  that  a*  is  a  minimizer  of  Problem  (7)  if  and  only  if  it  is  a  minimizer  of  Problem 
(8).  Moreover,  every  minimizer  can  be  written  as  a*  =  A^  (v  —  b)  +  w,  where  v  =  v  }e£^b/n)  ^  the 

pseudoinverse  of  A,  w  satisfies  Aw  =  0  and  0<  a*  <  1. 


We  next  show  a  relation  between  min-max  partitions  and  min  partitions. 

Theorem  5  (Min-max  and  min  partitions)  Let  £*2  be  a  m*n  partition  solution  to  Problem  2.  Then,  at,  is  also 
a  solution  to  Problem  1,  that  is, 


H-AaJj  +  &||oo  =  min_||Aa-|-&||00. 

OL<.OL<.Ot 


In  order  to  prove  Theorem  5,  we  introduce  the  following  definitions  and  results.  Given  a  partition  V  =  {V\,  ■  ■  ■ , Vn } 
defined  by  a ,  the  maximal  graph  associated  with  a  is  £/max  =  (V™ax,f“ax),  where  V™ax  =  {vi  €  Vc  :  L,  = 
max{Li, . . . ,  Ln}},  and  £cmax  =  (Vcmax  x  Vcmax)  n  8. 

Lemma  6  (Maximal  graph)  Let  a*  be  a  min  partition  of  the  graph  Q  =  (V,  £),  and  let  Qmax  =  (V™1,  £l™ax)  be 
the  maximal  graph  associated  with  a*.  Then,  for  each  n,  €  V™ax,Vj  £  Vc\  V™ax  with  {vi,vj}  £  8C,  it  holds  a*j  = 
if  i  <  j,  and  aT  =  if  i  >  j. 


PROOF.  Let  a*  be  a  min  partition,  and  let  L*  =  Aa*  +  b.  By  definition,  Vcmax  is  the  set  of  vertices  v,  such 
that  L*  =  || U*  || 00 -  Let  Vi  £  Vcmax,  and  partition  its  neighbor  set  as  Af)n  =  A f}  U  A/"2,  where  Nf  =  A Qn  (~l  Vcmax 
and  Mf  =  \  Nf .  Suppose  by  contradiction  that  af  >  for  some  Vj  £  A f}  with  i  <  j.  Define  a  from  a*  by 

modifying  only  the  entry  =  a*t]  —  e,  with  e  £  K>o-  Let  L  =  Aa  +  b ,  and  let  e  be  such  that  aij  >  ai;j  and  Li  >  Lj. 

An  equivalent  condition  for  Li  >  Lj  is  Ci  —  C2  >  £ij(a*i  —  a*j  +  2e),  where  ci  =  L*  —  ci  =  L*  —  a*tIij,  and 

=  1  —  a,ji  for  i  >  j.  It  follows  that 

WWl  -  I|i|||  =  (L*r  -  L\  +  (L*)2  -  L2j  =  (d  +  a*jlijf  -  (d  +  aijiijf  +  (c2  +  a*^)2  (c2  +  a^f 

=  2e%(ci  -c2  +  £ij  (a*j  -  -  e))  >  2e%  >  0, 

which  contradicts  our  assumption  of  a*  being  a  min  partition.  We  conclude  that  a,?  =  al3.  The  case  of  a*j  =  al3  is 
treated  analogously,  and  the  theorem  follows.  □ 
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We  are  now  ready  to  prove  Theorem  5. 


PROOF.  Let  a*  be  a  min  partition.  Recall  from  Lemma  6  that  there  exists  a  set  of  cameras  Vcmax  such  that  (i) 
L*  =  j|L*||oo  for  all  vt  G  Vcmax,  and  (ii)  a*j  =  atj  (resp.  a*j  =  5^)  if  i  <  j  (resp.  i  >  j)  for  all  {viyVj}  G  £c  with 
Vi  G  Vcmax  and  Vj  G  Vc  \  Vcmax.  Notice  also  that,  because  of  the  cameras  constraints  and  property  (ii),  it  holds 

min_  V  (Aa  +  6)i>|Vcm“|  Halloo. 

CK<  CK<  OL  Z - ' 


Let  a  be  such  that  ||L||oo  <  ||L*||oo-  Then  Li  <  L*  for  all  Vi  G  Vcmax-  It  follows  that 

(Ad  +  6)i<|Vcmax|||L*||00, 

fievcmax 


which  contradicts  our  hypothesis. 


□ 


Distributed  algorithms  to  compute  continuous  graphs  partitions  are  presented  in  Section  6.  In  the  next  section 
we  discuss  the  relationship  between  continuous  graph  partitions  and  the  design  of  trajectories  for  camera  network 
surveillance. 


3  Setup  for  Camera  Surveillance 

In  this  section  we  describe  our  setup  and  we  introduce  some  concepts  which  will  be  extensively  used  to  state  our 
results. 

3.1  Problem  Setup 

We  consider  the  problem  of  surveilling  an  environment  of  interest  by  means  of  a  camera  network.  We  represent 
the  environment  with  an  undirected  weighted  roadmap  Q  =  (V,£),  where  V  and  £  denote  the  vertex  and  the 
edge  sets,  respectively  [15].  In  particular,  each  vertex  Vi  G  V  corresponds  to  a  location  in  the  environment,  and 
{vi,Vj}  G  £  if  and  only  if  the  segment  [ty ,  Vj\  joining  vertices  ty  and  Vj  belongs  to  the  environment  (vertices  Vi 
and  Vj  are  within  line  of  sight).  Finally,  the  weight  of  the  edge  {vi,Vj}  G  £  equals  the  length  lij  =  ||i>j  —  Vj\\2,  and 
P1™  =  ma x{£ij  :  {viyVj}  G  £}. 

Cameras  are  placed  at  the  locations  Vc  C  V.  Define  the  set  of  neighboring  cameras  Mi  as  in  (1).  The  concept  of 
neighboring  cameras  will  be  exploited  in  Section  6  to  design  distributed  algorithms  for  the  cameras. 

Let  Xi(t )  denote  the  position  at  time  t  of  the  f.o.v.  of  the  *-th  camera.  We  assume  that  each  camera  has  a  limited 
visibility  range  along  each  adjacent  edge.  In  particular, 

(Al)  at  all  times  t ,  the  i-th  f.o.v.  is  a  point  along  the  segment  [ Vi ,  Vj }  for  some  Vj  G  Mp, 

(A2)  the  speed  of  the  i-th  f.o.v.  belongs  to  the  set  {0,1},  that  is,  the  f.o.v.  of  camera  Cj  either  is  stationary  at  some 
point  or  it  moves  at  maximum  (unitary)  speed; 

(A3)  for  each  Vj  G  A/),  a  point  Vij  G  [vi,Vj]  is  given  such  that  Xi(t)  G  [vi,Vij]  at  all  times  t; 

(A4)  the  cameras  locations  set  Vc  satisfies 


U  {{vi’vi}  :  vj  G  V}  =£, 

ViGVc 


so  that  the  roadmap  Q  is  jointly  visible  by  the  cameras. 
Our  setup  is  illustrated  in  Fig.  1. 
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Fig.  1.  This  figure  shows  an  environment  to  be  surveilled  by  a  camera  network.  The  environment  is  represented  by  a  roadmap 
Q  =  (V,£)  with  V  =  {ui, . . . ,  V14}.  Edges  £  are  denoted  with  solid  black  lines.  Cameras  are  installed  at  the  locations 
Vc  =  {vi, . . .  ,v?}.  White  rectangles  along  the  edges  represent  cameras  visibility  constraints,  and  the  parameters  a.  define  a 
continuous  partition  of  Q.  Finally,  the  DF-Trajectory  associated  with  the  partition  given  by  a  is  identified  by  the  closed  paths 
around  the  cameras. 

3.2  Cameras  trajectory 

A  cameras  trajectory  is  a  set  of  n  continuous  functions  X  =  {xi, . . . ,  xn},  where  Xj  :  K>0  — >•  £  describes  the  position 
of  the  z-th  f.o.v.  along  the  roadmap  Q .  We  focus  on  periodic  cameras  trajectories,  for  which  there  exists  a  finite  time 
T  £  R>o  satisfying  X(t  +  T)  =  X(t)  for  all  t  G  K>o-  Define  the  image  of  the  z-th  camera  as  the  set  of  points  visited 
by  the  z-th  f.o.v.  in  any  period  of  length  T,  i.e., 

Im(xi)  =  Ute[0jT] 

and  the  cameras  image  set  as  Ix  =  {Im(ii), . . . ,  Im(xra)}. 


Trajectory  1  DF-Trajectory  for  z-th  camera 
Input:  Parameters  cti j  and  set  of  neighbors  A/i; 

Set  Si(t j  —  (tVij  “1“  ( OCij  £%j  t )  'C'l )  j  O-'/J  £ij  ,  for  I,  £  [0,  OCij  £ij  ] , 

Set  to  =  0; 

for  Vj  £  Mi  do 

Set  Vij  as  in  Eq.  (3); 

Xi(t)  =  Si(t  -  t0),  for  t  £  [t0,  t0  +  atijtij]; 

Xi(t )  =  Si(2aij£ij  -  (t  -  t0)),  for  t  G  [to  +  oiij£ij,to  +  2 aij£ij]\ 
to  =  to  +  2  oiij£ijm, 

end  for 


We  now  define  a  particular  cameras  trajectory  associated  with  a  continuous  roadmap  partition.  The  optimality 
properties  of  this  trajectory  will  be  shown  in  the  subsequent  sections.  Let  o:df  define  the  continuous  partition  Vdf 
as  in  (2).  The  DF-Trajectory  Xdf  with  image  set  Vdl  is  obtained  by  letting  each  camera  sweep  its  subroadmap  in  a 
depth-first  order  [10],  and  it  is  formally  described  in  Trajectory  1.  See  Fig.  1  for  a  graphical  illustration. 

3.3  Performance  criteria 

In  this  work  we  design  cameras  trajectories  to  detect  intruders  along  the  roadmap.  We  consider  both  static,  and 
dynamic  intruders.  The  trajectory  of  an  intruder  is  a  continuous  function  p  :  R>o  — >  £.  Let  Ltd  be  the  set  of  all 
possible  intruder  trajectories,  and  let  IIS  the  set  of  static  intruder  trajectories  (pit)  =  po  for  all  t  >  to  and  for  some 
Po  €  £)■ 


An  intruder  is  detected  as  soon  as  its  position  coincides  with  the  f.o.v.  of  a  camera.  We  define  the  worst-case 
detection  time  of  a  cameras  trajectory  as  the  longest  time  for  the  detection  of  an  intruder.  In  particular,  for  an 
intruder  appearing  at  time  to  and  moving  with  trajectory  p,  and  a  cameras  trajectory  X,  let 

t*(t0,p,X)  =  min{{f  -  t0  :  t  >  t0,  p(t)  £  X(t)}  U  {00}}  . 
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We  define  the  static  worst-case  detection  time  as 


WDTS(X):=  sup  t*(t0,p,X),  (9) 

pens,  t0e[o,T] 

and  the  dynamic  worst-case  detection  time  as 

WDTd(X):=  sup  t*(t0,p,X).  (10) 

pend,  t0e[o,T] 


For  the  ease  of  notation  we  define 


WDT* 


inf  WDTS(X),  WDT*d  =  inf  WDTd(X). 

X  G  X 


We  conclude  this  section  by  observing  that  WDTS(X)  <  WDTd(X)  for  any  cameras  trajectory  X ,  and  that  for  any 
periodic  cameras  trajectory  X,  WDTS(X)  <  oo  if  and  only  if  the  entire  roadmap  is  persistently  surveilled  by  the 
cameras,  i.e. ,  £  C  Im(7fj.  Necessary  and  sufficient  conditions  for  a  trajectory  to  have  finite  dynamic  detection  time 
are  discussed  in  Section  5. 


4  Camera  Trajectory  for  Static  Intruders 

This  section  contains  our  results  for  the  detection  of  static  intruders. 

4-1  Main  results  for  static  intruders 

Consider  a  roadmap  Q  =  (V,£)  with  cameras  locations  Vc-  Let  V*  =  {Pf, . . . ,  V*\  be  a  continuous  partition  of  Q  of 
cardinality  |VC|  =  n  with  smallest  dimension,  that  is 

max  L*  =  min  max  Li,  (11) 

where  V  =  { V\ , . . . ,  Vn}  is  a  continuous  partition  of  Q.  Let  X*  be  the  DF-Trajectory  associated  with  the  partition 
V*.  Recall  that  the  roadmap  Q  with  cameras  locations  Vc  is  a  tree  if  the  induced  graph  Q(VC)  contains  no  cycles, 
and  it  is  a  ring  if  G(VC)  consists  of  a  single  cycle  [10]. 

Theorem  7  (Static  worst-case  detection  for  DF-Trajectory)  Consider  a  roadmap  Q  =  (V,£)  with  cameras 
locations  Vc  and  |VC|  =  n.  Let  X  be  the  DF-Trajectory  associated  with  a  continuous  partition  V*  of  Q  with  smallest 
dimension.  Then, 

(i)  WDTS(X*)  =  2 C(V*)  ,  and 

(ii)  WDTS(X*)  <  2 WDT*. 

Moreover,  if  Q  is  a  tree  or  a  ring,  then  WDTS(X*)  =  WDT*. 


In  Theorem  7  we  show  that  cameras  trajectories  designed  from  a  continuous  roadmap  partition  achieve  detection 
performance  within  a  constant  factor  of  optimal.  Since  optimality  is  guaranteed  for  tree  and  ring  roadmaps,  we  state 
the  following  conjecture. 

Conjecture  8  (Optimality  for  cyclic  roadmaps)  Motivated  by  our  results  in  Theorem  7,  we  conjecture  that 
WDTS(X*)  =  WDT*  also  for  cyclic  roadmaps. 
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4-2  Proof  of  Theorem  1 


In  this  section  we  derive  a  proof  of  Theorem  7.  We  start  by  introducing  the  necessary  notation  and  some  preliminary 
results.  We  define  a  relative  order  among  the  cameras  as  follows.  Let  {ci,...,cn}  be  the  set  of  cameras.  For  the 
neighboring  cameras  Ci  and  Cj  let  the  time  £  be  such  that  the  i-th  f.o.v.  and  the  j-th  f.o.v.  lie  on  the  edge  {vi,vj}. 
Then  we  say  that  Ci  <  Cj  if  ||x,;(£)  —  i>j||2  <  \\xj(t)  —  u»||2-  If  %i{t)  and  xj(t)  He  on  different  edges  at  time  £,  our 
convention  is  Xj(£)  <  Xj(t)  for  i  <  j.  A  camera  trajectory  is  order-invariant  if  the  relative  order  of  the  cameras  is 
preserved  over  time. 

Theorem  9  (Order-invariant  cameras  trajectory)  Given  a  roadmap  Q_=  (V,£)  with  cameras  locations  Vc  and 
a  cameras  trajectory  X ,  there  exists  an  order-invariant  cameras  trajectory  X  with  WDTS(A)  =  WDTS(A). 

PROOF.  Let  Ci  and  Cj  be  adjacent  cameras,  and  assume  the  existence  of  t  >  0  such  that  Xi(t)  =  Xj(t).  Define 
t^j  =  min{£  >  0  :  Xj(£)  =  Xj(t)}  and,  recursively,  ££•  =  min{£  >  :  Xi(t)  =  Xj(t)},n  G  N.  An  order-invariant 

trajectory  can  be  derived  from  X  permuting  the  cameras  labels  as  follows  (At  =  0,1,.. .): 

Xi(t)  =  Xj(t)  and  Xj(t)  =  Xi(t)  if  tfk  <  £  <  tfk+1, 

Xi(t)  =  Xi(t)  and  Xj(t)  =  Xj(t )  otherwise. 

Each  point  along  the  roadmap  is  visited  at  the  same  times  in  X  and  X,  thus  WDTS(X)  =  WDTS(X).  □ 


In  general,  the  images  of  neighboring  cameras  may  overlap.  A  camera  trajectory  is  called  non-overlapping  if  for 
every  pair  Cj  and  Cj,  it  holds  Int(Im(xj))  n  Int(Im(xj))  =  0,  where  Int  denotes  the  interior  of  a  set. 

Theorem  10  (Non-overlapping  cameras  trajectory  for  tree  and  ring  roadmaps)  Given  a  tree  (resp.  ring) 
roadmap  Q  =  (V,£)  with  cameras  locations  Vc  and  a  cameras  trajectory  X,  there  exists  an  order-invariant  and 
non-overlapping  cameras  trajectory  X  with  WDTS(A)  <  WDTS(X). 


PROOF.  Without  affecting  generality,  we  assume  that  the  trajectory  X  is  order-invariant  (cf.  Theorem  9).  We 
start  by  considering  tree  roadmaps,  and  we  define  the  trajectory  X  from  X  as  follows. 

Let  a i  =  |A/im|  +  1,  af)  =  |A/)ln|  for  i  =  2 ,n  and  let  Vf  =  0  for  every  i  =  1, . . .  ,n.  Iteratively  perform  the 
following  operations: 

Sk  =  {Vi  G  Vc  :  at1  =  1}, 

Vk  =  Cl(Im(xj)  \  (u^g^mp*"1)),  for  each  Vi  G  Sk, 
ak  =  ak_1  —  1,  for  each  Vj  G  Af™  U  {u,}  , 

where  Cl(-)  denotes  the  closure  of  a  set,  and  k  =  1,2,....  After  a  finite  number  kf  of  iterations,  the  set  Vkt  = 
{Vki, . . . ,  Vkt}  is  a  continuous  partition  of  Q.  Finally,  define  the  trajectory  X  as  the  DF-Trajectory  on  the  partition 
Vki .  In  the  interest  of  space,  we  do  not  prove  the  convergence  of  the  above  procedure,  and  we  provide  instead  an 
illustration  of  the  final  partition  in  Fig.  2. 

We  now  show  that  WDTS(A)  <  WDTS(A).  Consider  camera  c*,  let  |Mm|  =  rij,  and  let  d(fPkt)  =  {vn, . . .  ,ViUi}  be 
the  boundary  points  of  Vkl .  By  construction,  there  is  only  one  boundary  point,  say  u,;2 ,  such  that  some  points  in 
the  interior  of  [uj,  v^]  may  be  visited  by  a  camera  adjacent  to  Ci  with  reference  to  trajectory  X  (cf.  Fig.  2).  Notice 
that  every  boundary  point  is  visited  by  Cj.  Without  affecting  generality,  and  possibly  after  relabeling  the  boundary 
points,  let  t,\  be  such  that  Xi{t\)  =  Vn,  and  £2  such  that  ^(£2)  =  Va  and  no  other  boundary  point  is  visited  in 
Int([£i, £2]).  Notice  that  £2  —  £1  >  ||fj  —  Dilb  +  ||^i  —  D2II2)  and  that  every  boundary  point  must  be  visited  in  the 
interval  [£2,  £i+WDTs(A)].  It  follows  WDTS(X)  >  2  ll^i —  Dill 2,  and  this  must  hold  for  each  camera  c*.  Finally 
notice  that  WDTS(A)  =  maxi6p 2  J2j=i  ll1'*  —  vij lb)  so  that  the  statement  follows.  The  case  of  a  ring  roadmap 
can  be  treated  analogously.  □ 


Fig. 
We  now 


2.  This  figure  shows  a  tree  roadmap,  and  it  illustrates  the  partitioning  procedure  in  the  proof  of  Theorem  10. 
prove  Theorem  7. 


PROOF.  Statement  (i)  follows  from  the  definition  of  DF-Trajectory,  because  each  camera  sweeps  its  assigned 
subroadmap  at  maximum  speed  along  a  depth-first  tour. 


Regarding  statement  (ii),  consider  a  min  partition  a,  and  let  C/max  =  (Vcmax,  fcmax)  be  its  associated  maximal  graph 
(see  Lemma  6).  Define  Length(£max)  =  (Cfu;  t>j}e£c  llw*  ~ (i) (ii) * *  vj\\2’  anc^  notice  that 


WDT*|g„iax  > 


Length(t/max) 

|Vcmax|  ’ 


where  WDT*|e„ax  denotes  the  smallest  static  worst-case  detection  time  for  C/max.  Indeed,  since  =  Lj  for  each 
Vi,  Vj  £  Vmax,  each  camera  needs  to  sweep  (at  unitary  speed)  an  image  of  length  Lc"y  „aX| — -  for  C?max  to  be  covered. 
Moreover,  due  to  Lemma  6,  cameras  outside  C7max  cannot  visit  any  point  in  the  interior  of  t/max.  It  follows  that 


WDT*  >  WDT*|6 


Finally,  since  WDTS(A'*)  =  2L™y1)i£e| — -,  we  conclude  that  WDTS(A*)  <  2  WDT*. 

Consider  a  tree  (resp.  ring)  roadmap.  Due  to  Theorem  10,  there  exists  an  order-invariant  and  non-overlapping 
trajectory  X  with  WDTS(A')  =  WDT*.  To  conclude  the  proof,  we  have  WDTS(A'*)  <  WDTS(X)  since  Tx  is  a 
continuous  partition.  □ 


5  Cameras  trajectories  for  dynamic  intruders 

In  this  section  we  design  cameras  trajectories  for  the  detection  of  dynamic  intruders.  We  start  by  characterizing 
a  necessary  and  sufficient  condition  on  the  cameras  locations  for  the  existence  of  trajectories  with  finite  dynamic 
detection  time. 

Theorem  11  (Existence  of  trajectories  with  finite  dynamic  detection  time)  Given  a  roadmap  Q  =  (V,£) 
with  cameras  locations  Vc,  the  following  statements  are  equivalent: 

(i)  There  exists  a  cameras  trajectory  X  satisfying  WDTd(A)  <  oo; 

(ii)  For  every  Vi  £  Vc  with  |A/)|  >  3,  there  exists  Vj  £  Mfn  with  ai:j  =0  if  i  <  j  and  aij  =  1  if  i  >  j. 

The  following  result  is  useful  to  prove  Theorem  11. 

Lemma  12  (Finite  dynamic  detection  time  for  single  camera)  Given  a  roadmap  Q  =  (V,£)  with  cameras 

locations  Vc  =  {fi},  there  exists  a  cameras  trajectory  X  with  WDTd(A)  <  oo  if  and  only  if  |A/”i|  <  2. 
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PROOF.  To  show  sufficiency,  let  |A/i  |  <  2,  and  note  that  Q (Vc)  is  a  chain.  Let  X\  be  such  that  camera  c\  continuously 
sweeps  the  chain,  and  note  that  WDTd(A)  <  oo. 

To  show  necessity  of  the  statement,  notice  that  if  |A/i|  >  2,  an  intruder  may  choose  its  trajectory  p  so  that  p{t)  ^  v\ 
and  p(t  +  e)  €  [v\,Vj]  whenever  X\(t)  =  v\  and  x\ (t  +  e)  G  [ui,  Ufe],  for  e  G  K>o,  Vj,  Vk  G  A/"°ut,  and  j  ^  k.  □ 


We  are  now  ready  to  prove  Theorem  11. 


PROOF.  In  order  to  show  that  (ii)  is  a  necessary  condition  for  (i),  suppose  that  (ii)  does  not  hold.  Then  camera 
Cj  needs  to  surveil  a  subroadmap  in  which  \J\ft  >  3.  The  statement  follows  from  Lemma  12. 

We  now  show  that  (ii)  is  also  a  sufficient  condition  for  (i)  by  proposing  a  procedure  to  clear  every  subroadmap  from 
intruders  appearing  at  time  0.  By  periodically  repeating  this  procedure,  intruders  appearing  at  different  times  are 
also  detected.  Notice  that  intruders  appearing  along  the  edge  [u*,  Vj\  can  be  detected  by  moving  the  cameras  c,  and  Cj 
towards  each  other  in  a  way  that  Xi(t)  =  Xj(t)  for  some  time  t.  Starting  from  iq,  if  |A/i|  <  2,  then  sweep  the  adjacent 
edges  by  synchronizing  the  motion  of  Vi  and  its  neighboring  cameras  in  any  order.  If  |A/i  |  >  2,  suppose  {iq,  U2}  G  Vc, 
then  let  a12  =  0  from  condition  (ii).  Let  camera  C2  sweep  the  entire  edge  {/tq , V2}  and  stop  at  iq,  and  let  the  other 
neighboring  cameras  Cj  G  stop  at  their  vertices  Vj.  Then  camera  c\  sequentially  sweeps  its  assigned  adjacent 
segments  [tq,iq j],  with  {vi,Vj}  G  £c,j  7^  2  by  synchronizing  its  motion  and  those  of  the  neighboring  cameras,  while 
keeping  the  f.o.v.  of  C2  at  vertex  tq.  Notice  that  any  intruder  appearing  in  an  edge  adjacent  to  v\  at  time  0  is  detected 
by  this  procedure.  Then  let  Ci  and  its  neighboring  cameras  return  to  their  vertices.  Finally  iterate  the  procedure  for 
subsequent  cameras  in  increasing  order,  and  repeat  over  time.  □ 


In  the  remainder  of  this  section  we  design  cameras  trajectories  for  the  special  cases  of  ring  and  tree  roadmaps.  We 
refer  the  reader  to  [21]  for  a  solution  to  the  case  of  chain  roadmaps  (cf.  Equal- Waiting  trajectory). 

5.1  Ring  roadmap 

Consider  a  ring  roadmap  Q  =  (V,  £)  with  cameras  locations  Vc,  and  visibility  constraints  a,  a.  Let  V*  =  {V*, . . . ,  V*} 
be  a  min-max  partition  of  Q,  and  let  L  =  C(V*)  =  max  {Li, . . . ,  Ln},  where  L j  is  defined  in  (4).  Notice  that 
\J\f-n\  =  2  for  all  Vi  G  Vc,  and  that  each  subroadmap  V*  can  be  written  as  a  segment,  parametrized  by  s,  :  [0,  L,]  — ► 
,  such  that  s,;(t)  —  {I’Vi^i^- 1  T  (L^  f)u^_r,i)/Lj. 

Trajectory  2  Ring-Sync- Trajectory  for  camera  Cj 

Camera  a _surveils  V*  =  w»]  U  [vi,Vi,i+ 1]; 

Set  In  =  [L,  2 L\,Tn  =  L,  if  n  is  odd; 

In  =  0,  Tn  =  0,  if  n  even; 
if  i  is  even  then 

Xi(t)  =  for  t  G  [0,  L  -  Li];  _ 

Xi(t)  =  Si(t  —  (L  -  Li)),  for  t  G  [L  —  Li,L\, 

Xi(t)  =  Vj.i+i,  for  t  G  In  U  [L  +  Tni2L  -  Li  +  Tn ]; 

Xi(t)  ~  Si(2L  -(-  Tn  —  t ),  for  t  G  [2 L  —  Li  -(-  Tn,  2 L  -t-  Tn] ; 
else  if  i  is  odd  then 

Xi(t)  =  vi:i+ 1,  for  t  G  [0,  L  —  Li]; 

Xi(t)  =  Si(L  -  t),  for  t  G  [L  —  Li,L];_ 

Xi(t)  =  Vi— i.i,  for  t  G  In  U  [L  +  Tn,  2 L  —  Li  -f  Tn]; 

Xi (t)  —  Si(t  —  (2 L  —  Li  -f-  Tn)),  for  t  G  [2 L  —  Li  +  Tn ,  2 L  +  Tn] ; 
if  i  =  n  then 

Xi(t)  =  Vi,i+ 1,  for  t  G  [0,  2 L  —  Li]; 

Xi(t)  =  Si  (2L  —  t) ,  for  t  G  [2L  —  Li ,  2L] ; 

Xi(t )  =  Vi-i,i,  for  t  G  [2 L,  3 L  —  Li];_ 

Xi(t)  =  Si(t  —  (3 L  -  Li)),  for  t  G  [3L  —  Li,3L]; 

end  if 
end  if 


We  propose  the  Ring-Sync-Trajectory  in  Algorithm  2;  see  Fig.  3  for  a  graphical  illustration.  We  next  provide  an 
informal  description  of  the  Ring-Sync-Trajectory  within  its  period  for  the  case  of  an  even  number  of  cameras: 
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Fig.  3.  Consider  the  ring  roadmap  Q  =  (V,  £)  with  cameras  locations  Vc  and  n  =  5.  Its  associated  optimal  partition  V  depicted 
in  Fig.  3(a).  The  corresponding  3 L-  periodic  Ring-Sync- Trajectory  described  in  Trajectory  2  is  shown  in  Fig.  3(b). 

(i)  the  f.o.v.  of  camera  Cj  is  set  at  its  right  (left)  boundary  point  if  i  is  odd  (even),  so  that  all  the  cameras  are 
pairwise  synchronized; 

(ii)  each  camera  sweeps  its  subroadmap  at  maximum  speed,  and 

(iii)  each  camera  stops  for  L  ~  Li  at  each  boundary  point. 

If  the  number  of  cameras  is  odd,  the  algorithm  is  more  involved  and  it  is  formally  described  in  Trajectory  2.  Notice 
that,  (i)  the  i-th  camera  surveils  only  V* ,  (ii)  the  trajectory  of  each  camera  is  3L-periodic  if  n  odd,  and  2L-periodic 
if  n  even,  and  (iii)  the  f.o.v. s  of  each  pair  of  neighboring  cameras  coincide  at  some  times  within  the  period. 


Theorem  13  (Dynamic  worst-case  detection  for  Ring-Sync- Trajectory)  Given  a  ring  roadmap  Q  =  (V,£) 
with  cameras  locations  Vc,  let  Xs  be  the  Ring-Sync-Trajectory  in  Algorithm  2.  Then 


(i)  if  n  is  even, 

(ii)  if  n  is  odd, 


WDTd(Xs)  =  WDTS(XS)  =  WDTd; 
WDTd(Xs)  =  WDTS(XS)  <  ^WDTd. 


PROOF.  In  order  to  prove  statement  (i),  consider  a  ring  roadmap  with  n  even,  and  compute  the  optimal  partition 
V*  with  dimension  L  (see  Section  6).  Once  we  synchronize  the  cameras  according  to  Trajectory  3,  the  minimum 
dynamic  worst-case  detection  time  is  achieved.  In  order  to  prove  (ii),  notice  that  the  Ring-Sync- Trajectory  Xs  is 
3L-periodic  when  n  is  odd,  and  that  each  camera  surveils  the  corresponding  subroadmap  V*  ■  The  thesis  is  proved 
observing  that  WDTd  >  WDT*  =  2 L.  □ 


5. 2  Tree  roadmap 

Consider  a  tree  roadmap  Q  =  (V,£ )  with  cameras  locations  Vc  =  V  \  {tq},  where  node  iq  is  labeled  as  root.  Assume 
that  each  camera  can  entirely  surveil  each  adjacent  edge,  that  is,  a  =  0  and  5=1.  We  now  design  a  cameras 
trajectory  for  dynamic  intruders,  and  we  show  that  the  performance  of  our  trajectory  is  within  a  constant  factor  of 
optimality.  We  start  by  recalling  some  definitions  [10].  Then,  vertex  Vi  £  Vc  is  a  parent  of  vertex  Vj  £  Vc  ( Vi  is  a 
child  of  Vj)  if  {i’i,Vj}  £  £c  and  iq  lies  on  the  unique  shortest  path  from  Vj  to  tq.  Let  vf  denote  the  parent  of  iq,  and 
let  £?  denote  the  length  of  the  segment  [iq,  i>?].  Recall  that  £max  is  the  length  of  the  longest  edge.  Let  us  parametrize 
the  segment  [iq,  i;f]  as  Si  :  [0, £?]  — ►  [vi,  uf],  Si(t)  =  ( tvf  +  —  t)vi)/£?.  We  divide  the  vertices  into  two  groups 

according  to  the  parity  of  their  distance  to  the  root.  In  particular,  we  define  dist-r(uj)  as  the  number  of  edges  in  the 
shortest  path  from  tq  to  the  root  v\. 

We  propose  the  Tree-Sync-Trajectory  in  Trajectory  3  (see  Fig.  4  for  an  example).  An  informal  description  of  the 
Tree-Sync- Trajectory  follows: 
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Fig.  4.  Fig.  4(a)  show  a  tree  roadmap  where  v\  is  labeled  as  root  vertex.  In  Fig.  4(b)  we  report  the  Tree-Sync- Trajectory 
described  in  Trajectory  3.  Note  that  cameras  are  synchronized,  that  is,  for  each  pair  of  adjacent  cameras  {i>i,u,}  £  £c,  there 
exists  t  £  [0,  2fmax]  such  that  Xi(t)  =  Xj(t). 


Trajectory  3  Tree-Sync-Trajectory  for  camera  c* 

Camera  Ci  surveils  only  the  segment  [Vi,  Uip]; 
if  dist-r(ui)  is  odd  then 

Xi(l)  '=  for  t  £  [0,£max-£?]; 

Xi(t)  =  Si(fmax  -  f),  for  t  £  [fmax  -  fp,fmax]; 

Xi(t)  =  Vi,  for  t  £  [£max,  2fmax  -  fp]; 
aJi(i)  =  Si(t  -  (2tmax  -  fp)),  for  *  €  [2£max  -  fp,  2£max]; 
else  if  dist-r(ui)  is  even  then 

Xi(t)=Vi,  for  t  £  [0,  £max  —  ff]; 

a*(t)  =  Si(t  -  (fmax  -  ff)),  for  t  £  [fmax  -  £f,  fmax]; 

aJi(t)  =  Uip,  for  f  £  [£max,  2fmax  -  £p]; 

®i(i)  =  Si( 2fmax  -  t),  for  t  £  [2fmax  -  tp,  2f max] ; 
end  if 


(i)  the  f.o.v.  of  camera  c,  is  set  at  the  vertex  vf  (resp.  Vi)  if  dist-r(vj)  is  odd  (resp.  even), 

(ii)  camera  Ci  sweeps  the  segment  [vi,  vf],  and 

(iii)  camera  Ci  stops  for  £max  —  at  each  boundary  point. 

Notice  that,  (i)  camera  Cj  surveils  only  the  segment  [uj,  vf],  (ii)  the  trajectory  of  each  camera  is  2£max-periodic,  (iii) 
cameras  are  synchronized ,  that  is  Xi(tk )  =  xf  (t*,)  for  tk  =  fcfmax,  with  k  £  N  even  (resp.  odd)  if  dist-r(vj)  is  even 
(resp.  odd).  Our  Tree- Sync- Trajectory  in  Trajectory  3  extends  the  concept  of  Equal- waiting  trajectory  to  the  case 
of  tree  roadmap  [21]. 

Theorem  14  (Dynamic  worst-case  detection  for  Sync- Trajectory)  Given  a  tree  roadmap  Q  =  (V,£)  with 
cameras  locations  Vc  =  V  \  {vi}  and  at  =  0,  a  =  1,  let  Xs  be  the  Tree-Sync-Trajectory  in  Trajectory  3.  Then, 

WDTd(Xs)  =  WDTS(X(i) * * * S)  <  2WDTd. 


PROOF.  Consider  the  Tree-Sync-Trajectory  Xs.  Notice  that  Xs  is  2fmax-periodic,  that  each  camera  surveils  a  single 
edge,  and  that  each  edge  is  surveilled  by  a  different  camera.  Let  an  intruder  appear  at  time  to  along  [v*,  v,;p],  with 
|| v»  —  p(fo)||2  <  || Vi  —  Xi(to)||2-  Then,  such  intruder  is  confined  in  [vj,  Vjp]u(uUj.  [v*,  Vj])  to  avoid  detection,  where  Vj  is  a 
child  of  Vi.  Since  vt  and  all  its  children  synchronize  at  most  every  2fmax,  we  have  WDTd(Xs)  =  WDTS(XS)  =  2fmax. 
Notice  that  the  case  ||vj  —  p(0)||2  >  ||vj  —  Xj(to)||2  yields  the  same  conclusion.  On  the  other  hand,  since  each  edge 
can  be  surveilled  by  at  most  two  cameras,  it  holds  WDTd  >  WDT*  >  ^max;  which  concludes  the  proof.  □ 


6  Distributed  Partitioning  Algorithms 

In  what  follows  we  design  three  distributed  algorithms  for  the  continuous  min-max  partitioning  problem.  Given  an 
optimal  partition,  cameras  organize  along  a  DF-Trajectory  as  in  Trajectory  1. 
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Algorithm  4  Synchronous  Gradient  Partitioning 

for  Vi  G  Vc  do 

Camera  d  receives  Sj  from  Cj,  for  all  Vj  €  A/j‘n; 
Set  ajf1  <-  afj  -  -  L]); 

if  ajt1  <  atj  then  a*tl  =  aij; 
else  if  aft1  >  ouj  then  a-t1  =  aii; 

end  if 

Camera  a  transmits  S’*"1"1  to  Cj,  for  all  Vj  £  A/*". 

end  for 


The  algorithms  we  present  rely  upon  different  cameras  communication  assumptions.  We  assume  each  camera  to  be 
equipped  with  a  wireless  sensor  device.  In  all  our  algorithms,  cameras  perform  the  following  operations:  (i)  receive 
parameters  from  (some)  neighboring  cameras,  (ii)  update  the  parameters  corresponding  to  (some)  adjacent  edges, 
and  (iii)  transmit  the  new  values  to  (some)  neighboring  cameras.  These  operations  are  detailed  in  the  next  sections. 
For  convenience,  let  S\  =  {aP  :  vj  £  A/-)11}  be  the  state  of  camera  C;  at  iteration  t  £  N.  Finally,  we  initialize 
aij  =  f°r  all  {vu  vj}  G  £c  with  i  <  j. 


6.1  Synchronous  Gradient  Partitioning  algorithm 

The  distributed  algorithm  presented  in  this  section  assumes  a  synchronous  mode  of  operation  of  the  cameras,  and 
it  is  inspired  by  the  classical  gradient  projection  method  [5].  In  particular,  every  camera  performs  operations  at 
uniform  time  instants.  The  f-th  iteration  of  this  algorithm  is  detailed  in  Algorithm  4. 

Theorem  15  (Synchronous  Gradient  Partitioning)  For  a  roadmap  Q  with  cameras  locations  Vc,  let  A  and 
b  be  as  in  (5)  and  (6),  respectively.  Let  0  <  e  <  ((|maxfnax)-1)  where  dmax  =  max  { |  A/-)11 1  :  Vi  6  Vc}-  Then,  the 
Synchronous  Gradient  Partitioning  algorithm  in  Algorithm  4  asymptotically  converges  to  agGD  =  limt_).0O  a*,  and 

min  ||  Aol  +  6||^,  =  ||Ac*gGD  +  6)1^,, 

a<a<a 


where  a  and  a  denote  the  cameras  constraints. 


PROOF.  Note  that  the  update  step  can  be  expressed  in  vector  form  as 

OLt+1  t—  ct4  —  £AT(Aat  +  b), 

and  that  ATAct  +  ATb  is  the  gradient  of  the  quadratic  function  |||Aa  +  6|||.  Therefore  the  Synchronous  Partitioning 
algorithm  coincides  with  the  gradient  projection  method  [5].  To  conclude  the  proof  note  that  the  gradient  of  -)||Aa:  + 
bill  is  Lipschitz-continuous  with  some  Lipschitz  constant  K  G  R>o-  Thus,  for  a  sufficiently  small  step  size  e,  precisely 
if  0  <  e  <  2/K,  the  convergence  of  cA  to  a  min  partition  is  guaranteed  by  [5,  Proposition  3.4].  In  order  to  compute 
the  stated  upper  bound  for  the  stepsize  e,  it  can  be  shown  that  K  <  2dmax(tmax)2.  Finally,  the  claimed  statement 
follows  from  Theorem  5.  See  [5]  for  further  details.  □ 


6.2  Asymmetric  Broadcast  Partitioning  algorithm 

The  distributed  algorithm  presented  in  this  section  assumes  an  asymmetric  broadcast  communication  protocol.  In 
particular,  at  each  iteration  only  one  camera  updates  its  state  by  using  local  information  from  its  neighboring 
cameras.  In  order  to  guarantee  the  convergence  of  the  algorithm,  we  assume  the  existence  of  a  finite  duration 
r  G  R>o  such  that,  for  all  t  G  R>o,  every  camera  in  Vc  is  selected  at  least  once  in  the  time  interval  [t,  t  +  r)  ( partial 
asynchronism  assumption).  The  f-th  iteration  of  this  algorithm  is  detailed  in  Algorithm  5. 

Theorem  16  (Asymmetric  Broadcast  Partitioning)  For  a  roadmap  Q  with  cameras  locations  Vc  and  cameras 
edges  £c,  let  A  and  b  be  as  in  (5)  and  (6),  respectively.  Let  r  be  the  partial  asynchronism  constant,  and  let  0  < 
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Algorithm  5  Asymmetric  Broadcast  Partitioning 
Camera  d  is  randomly  selected; 

Camera  d  receives  S!j  from  camera  Cj,  for  all  Vj  £  A /?“; 
for  Vj  £  A/)1"  do 

<4+i  <-  ah  -  — Lj)  i 

if  a-/1  <  a;j  then  a*/1  =  atj; 
else  if  a**1  >  aij  then  aV1  =  di-,; 

end  if 
end  for 

Camera  d  transmits  S(+1  to  camera  Cj,  for  all  Vj  £  A Qn. 


Algorithm  6  Symmetric  Gossip  Partitioning 

Neighboring  cameras  d  and  c,  are  randomly  selected; 
Camera  d  ( Cj )  receives  S)  (S’*)  from  Cj  (c;); 

L*  =  {Ll  +  L))/ 2; 

aij  =  {L  ~^2vkeA 

if  ajj-1  <  atj  then  a*/1  =  a^; 

else  if  a*+x  >  a ij  then  =  a ij; 

end  if 

Camera  d  transmits  S*+1  to  Ck,  for  all  vt  £  A Qn; 
Camera  Cj  transmits  S*+1  to  Ck,  for  all  Vk  £  A /”jn. 


£  <  (K(l  +  t  +  t|Sc|))  1 ,  where  K  £  R>o  is  the  Lipschitz  constant  of  a  — ►  AT(Aa  +  b).  Then,  the  Asymmetric 
Broadcast  Partitioning  algorithm  in  Algorithm  5  asymptotically  converges  to  c*ab  =  limt->oo  of .  Moreover , 

min_||Aa  +  6||^  =  ||Aa^B  +  6||^, 

ot<ot<cx 


where  a  and  a  denote  the  cameras  constraints. 


PROOF.  As  in  Theorem  15,  the  algorithm  update  follows  the  gradient  of  a  — »  |||  Aa  +  b Hi-  Because  of  the  partial 
asynchronism  assumption  and  the  fact  that  a*  £  [0, 1]£°  is  such  that  a  <  a*  <  a  for  all  t  £  N,  the  statement  follows 
from  [5,  Section  7,  Proposition  5.3]  and  Theorem  5.  Note  that  the  bound  for  the  stepsize  e  depends  on  the  Lipschitz 
constant  K,  the  time  horizon  r  and  the  number  of  edges  connecting  cameras  £c  \ .  See  [5]  for  further  details.  □ 


6.3  Symmetric  Gossip  partitioning  algorithm 

The  distributed  algorithm  presented  in  this  section  assumes  a  symmetric  gossip-type  communication  protocol.  In 
particular,  at  each  time  iteration  only  one  component  of  a  camera  state  is  updated,  and  only  two  adjacent  cameras 
are  involved  in  the  computation.  The  f-th  iteration  of  this  algorithm  is  detailed  in  Algorithm  6. 

Theorem  17  (Symmetric  Gossip  Partitioning)  For  a  roadmap  Q  with  cameras  locations  Vc,  let  A  and  b  be  as 
in  (5)  and  (6),  respectively.  Let  the  partial  asynchronism  assumption  hold.  Then,  the  Symmetric  Gossip  Partitioning 
algorithm  in  Algorithm  6  asymptotically  converges  to  ccgG  =  lim^oo  of .  Moreover, 

min_||Aa  +  6||^  =  ||  AagG  +  6||^, 

ot<ot<ot 


where  a  and  a  denote  the  cameras  constraints. 


PROOF.  Define  17(a)  =  E{v,,Vj}e£c  (Li  —  Lj)2  as  energy  storage  function.  The  convergence  of  the  algorithm  can 
be  retrieved  reasoning  along  the  lines  of  [1,  Theorem  IV.  1],  and  by  applying  Theorem  5.  □ 
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Fig.  5.  This  figure  shows  the  convergence  of  the  Synchronous  Gradient  Partitioning  (SGD,  blue  solid  line),  the  Asymmetric 
Broadcast  Partitioning  (AB,  black  dash-dot  line),  and  the  Symmetric  Gossip  Partitioning  algorithms  (SG,  green  dashed  line) 
towards  a  solution  of  the  continuous  min-max  partitioning  problem.  For  the  simulation  we  use  the  configuration  in  Fig.  2, 
with  a  =  0,  a  =  1  and  a0  =  0. 

To  conclude  this  section,  we  validate  our  distributed  algorithms  through  a  numerical  study.  The  tree  roadmap 
considered  for  the  simulations  is  in  Fig.  2.  Notice  that  \£c\  =  6,  n  =  |VC|  =  7,  and  the  number  of  locations  is  given 
by  |V|  =  14.  The  stepsizes  for  Algorithm  5  and  Algorithm  6  are  chosen  sharp,  up  to  a  constant  e  =  .01,  to  their 
upper  bounds  stated  in  Theorem  15  and  Theorem  16,  respectively.  The  results  of  our  simulation  study  are  in  Fig. 
5.  Notice  that  all  the  proposed  algorithms  converge  to  the  desired  value. 

7  Conclusion 

In  this  work  we  design  surveillance  trajectories  for  a  network  of  autonomous  cameras.  We  consider  both  static  and 
dynamic  intruders  in  the  environment  to  be  surveilled.  As  performance  criteria  we  consider  the  worst-case  detection 
time  of  intruders.  For  the  case  of  static  intruders,  we  derive  optimal  trajectories  for  ring  and  tree  roadmaps,  and 
constant-factor  suboptimal  trajectories  for  general  roadmaps.  For  the  case  of  dynamics  intruders,  we  derive  constant- 
factor  suboptimal  trajectories  for  chain,  ring,  and  tree  roadmaps.  As  a  complementary  result,  we  introduce  the 
continuous  partitioning  problem  of  a  weighted  graph,  and  we  propose  distributed  algorithms  for  its  solution. 
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